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Abstract 

The anharmonic force field of BeH 2 has been calculated near the basis set and n- 
particle space limits. The computed antisymmetric stretch frequencies of BeH 2 and BeD 2 
are in excellent agreement with recent high-resolution gas-phase measurements [P. F. 
Bemath et ah, Science 297, 1323 (2002)]. The agreement between theory and experiment 
for the other spectroscopic constants is also excellent, except for co 3 and X 33 for BeH 2 and 
G 22 for BeD 2 . It is concluded that further experimental work is needed in order to resolve 
these discrepancies. 


Introduction 

Recently, Bemath and coworkers reported the first high-resolution rovibrational 
band assignment for polyatomic molecules containing beryllium, BeH 2 and BeD 2 [1,2]. 
Beryllium vapor is known to be very toxic, which probably explains why there have been 
very few gas-phase experimental studies on beryllium-containing polyatomic molecules to 
date. Nonetheless, BeH 2 has been the target of several theoretical studies [3-5] over the 
years, partly because it has only four valence electrons and thus is an ideal test molecule 
for new electron correlation methods. BeH 2 has been the subject of very accurate 
theoretical studies, such as Ref. [3], where a quartic force field, equilibrium geometry, and 
fundamental vibrational frequencies were computed using the CCSD(T) method [6] 
(singles and doubles coupled-cluster theory with a perturbational estimate of connected 
triple excitations). High-level ab initio calculations have also shown that BeH 2 is 
thermodynamically stable [4] relative to dissociation to Be + H 2 (37.6 kcal/mol), but 
direct insertion of a beryllium atom into the H 2 bond has a barrier of 48.6 kcal/mol [5], 

The v 2 and v 3 bands of BeH 2 and BeD 2 were first observed [7] in an argon matrix 
infrared (IR) study of reaction products of H 2 with beryllium atoms generated by pulsed 
laser evaporation. 

Bemath and coworkers’ assignment for the antisymmetric Be-H stretch, v 3 , was 
aided by the aforementioned CCSD(T) study of the fundamental vibrational frequencies 
[3]. In this case, the theoretical calculations predicted a fundamental vibrational 
frequency of 2167.2 cm' 1 , which compares reasonably well with the experimental 
assignment of 2178.8 cm' 1 . (The latter, incidentally, implies a sizeable Ar matrix red shift 
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of 19.7 cm' 1 in Ref. [7].) However, the calculations in Ref. [3] included only the valence 
electrons in the CCSD(T) correlation procedure and only considered up through /-type 
functions in the one-particle basis set. There is ample evidence now to indicate that the 
majority of the remaining error in valence CCSD(T) calculations is due to the neglect of 
core correlation, and this is even more likely to be the case for beryllium compounds due 
to the well known near degeneracy between the Is and 2s shells (for example, see Refs. 
[8-12]). The primary purpose of this study is to present state-of-the art calculations of 
the quartic force field, equilibrium structure, and fundamental vibrational frequencies for 
BeH 2 in order to assess the level of accuracy that is obtainable for polyatomic beryllium 
compounds from purely ab initio calculations. This study therefore includes the effects of 
core correlation and the use of one-particle basis sets that go beyond f - type functions, as 
well as valence correlation that goes beyond CCSD(T) and scalar relativistic 
contributions, of which the latter two are expected to be smaller effects. 

The theoretical approach is detailed in the next section followed by results and 
discussion. The final section contains our conclusions. 


Theoretical Approach 

Most electronic structure calculations have been carried out using the MOLPRO 
2002.3 electronic structure system [13]. The CCSD(T) [6] electron correlation method 
was used for most calculations. For n-particle calibration purposes, CCSDT (coupled 
cluster with all single, double, and triple excitations [14]) and FCI (full configuration 
interaction) calculations were also carried out, the former using the Austin-Mainz version 
of ACES II [15], the latter by means of the Handy-Knowles determinantal FCI code [16] 
in MOLPRO. 

Standard Dunning cc-pVnZ (correlation consistent polarized valence n- tuple zeta 

[17] ) basis sets were used throughout on hydrogen. Some calculations in which only 
valence electrons were correlated employed preliminary cc-pVnZ basis sets for beryllium 

[18] ; for calculations in which all electrons were correlated, the CVnZ (core-valence n- 
tuple zeta) basis sets for beryllium were taken from a recent paper on core-valence 
correlation in the alkali and alkaline earth metals [19]. 

Scalar relativistic effects were considered by means of the Douglas-Kroll 
approximation [20,21] as implemented in MOLPRO 2002.3. 

Spectroscopic constants for the molecule were obtained using standard second- 
order rovibrational perturbation theory [22] as implemented in the SPECTRO [23,24] 
program. 

Results and Discussion 
A. Spectroscopic Constants for BeH 2 

Computed and observed anharmonic spectroscopic constants are presented in 
Table 1 . Computed geometries and harmonic frequencies at some additional levels of 
theory can be found in Table 2. 

Using the cc-pVTZ basis set (i.e. As2p2d\f on Be and 3s2pld on H), we were able 
to carry out valence FCI calculations, which are an exact solution of the nonrelativistic 
Schrodinger equation in the given one-particle basis set. As expected for a molecule like 


- 2 - 


BeH 2 which is dominated by a single reference configuration, the CCSD(T) data with the 
same basis set are in veiy close agreement with FCI. The chief differences are that 
CCSD(T) apparently overestimates the three fundamental frequencies by 1.5, 0.4, and 1.1 
cm' 1 , respectively, and underestimates the bond distance by 0.00015 A. The largest 
anharmonicity constant, X 13 , is underestimated (in absolute value) by 0.13 cm' 1 . Higher- 
order effects of connected triple excitations account for essentially all of the difference: 
the CCSDT and FCI data are nearly indistinguishable. 

As the basis set is enlarged to 5s4p3d2f on Be and 4s3p2d on H (i.e., cc-pVQZ 
without g functions on Be nor / functions on H), the differences between CCSD(T) and 
FCI actually decrease for the frequencies, although the tiny geometric difference remains 
constant. This confirms that CCSD(T) is a suitable level of theory for treating the 
system. 

As expected, scalar relativistic effects are nearly nonexistent in this very light 
system: the only changes seen upon introducing Douglas-Kroll corrections are a further 
shortening of the bond distance by 0.0001 A and a decrease in the bending frequency by 
0.1 cm' 1 . 

As could be expected for this system based on previous experience with BeH [9], 
inner-shell correlation effects are quite prominent. Comparing CCSD(T)/CVQZ force 
fields with and without Befls) correlation, we find that the latter shortens the bond by 
about 0.004 A and increases the fundamental frequencies by 8, 3, and 9 cm' 1 , 
respectively. In fact, the CCSD(T)/CVQZ spectroscopic constants are in excellent 
agreement with the available experimental data: v 3 is reproduced to within 0.5 cm' 1 and r e 
to four decimal places, while X] 3 and X 23 agree to 0.25 and 0.06 cm' 1 , respectively. The 
fairly small discrepancy in G 22 might be eliminated if more experimental data involving v 2 
were available. Agreement is equally satisfying for the centrifugal distortion and 1- 
doubling constants. Disagreements (with opposite signs) of 0.001 cm' 1 are seen in oq and 
a 2 . These could reflect either contributions of higher-order rovibrational coupling in the 
experimental data, nonadiabatic effects (which are somewhat significant for the BeH bond 
distance [9]), or a combination of both. 

In light of the generally excellent agreement between theory and experiment — 
particularly for v 3 — the fact that our computed co 3 exceeds the experimentally derived 
value by nearly 5 cm' 1 is somewhat puzzling. The phenomenon appears to be linked with 
a discrepancy of -2 cm' 1 (theory predicting the larger absolute value) for X 33 . It is 
extremely unlikely that further improvement of the theoretical treatment will cause 
substantial changes in the computed X 33 . Furthermore, no resonances exist and the 
molecule is rigid enough for a second-order rotation-vibration perturbation theory 
treatment to be adequate. We note that the experimental value — obtained from 
[G(00°2)-2G(00°l)]/2 — is an effective one that includes higher-order anharmonicity 
contributions, X 33 eff =X 33 +Y 133 /2+Y 233 +9Y 333 /2+0(Zy k |). Y 333 can be estimated from 
[G(00 0 3)-3G(00 0 2)+3G(00°l)]/6=Y 333 -i-0(Zjji { |). If we assume the third anharmonicities to 
be negligible, then Y 333 is found to be -0.107 cm-1, nearly identical to the best computed 
co e y e for BeH [9]. This suggests either unusually large Y ]33 and Y 233 constants, or a 
problem with X 33 itself. Further experimental work would certainly be useful for 
clarifying this issue. 
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At the CCSD(T)/CVTZ level, v 2 is calculated to have an infrared intensity of 
472.8x2 km/mol, compared to 274.2 km/mol for the antisymmetric stretch which was 
observed. This certainly suggests the v 2 band ought to be amenable to experimental study 
in the gas phase. 

As noted previously [9,10], the anharmonicity constants are rather more sensitive 
to the basis set than generally assumed. Their sensitivity to inner-shell correlation is 
rather minor, reflecting mostly the changes in the underlying harmonic frequencies. Table 
1 suggests, however, that the CCSD(T)/CVQZ values are largely converged with the 
computational level of theory, and that observed values for the missing constants are 
unlikely to deviate substantially from our computational predictions. 

Further basis set expansion to CV5Z leads to a shortening of the bond by 0.0003 
A (half of which cancels with inclusion of higher-order correlation effects), and minor 
changes in the computed frequencies that largely cancel with inclusion of higher-order 
triple excitation effects. 

B. Spectroscopic Constants for BeD 2 

Agreement between CCSD(T)/CVQZ and experiment for the v 3 fundamental of BeD 2 
bolsters our conclusions about the quality of the BeH 2 force field. Indeed, the difference 
here is only 0.2 cm' 1 versus 2.1 cm' 1 for BeH 2 . For BeD 2 , the main disagreement with 
experiment is for G 22 , which is computed nearly an order of magnitude larger than the 
experimental value. Again, this suggests that further experimental work on the bend is in 
order, especially considering the excellent agreement for v 3 , a direct observable. 


Conclusions 

A highly accurate anharmonic force field for BeH 2 has been computed, and is 
capable of reproducing the observed data with near- spectroscopic accuracy. The principal 
shortcoming of the earlier calculation [3] was neglect of inner-shell correlation: the effect 
of higher-order connected triple excitations is an order of magnitude less important, and 
the effects of connected quadruple excitations and of relativity can be neglected outright. 
Some issues with the experimentally derived harmonic antisymmetric stretching 
frequency and the associated anharmonicity for BeH 2 suggest further experimental work 
may be in order, as does the 1-type resonance constant G 22 for BeD 2 . The bending 
frequency is predicted to be very intense and should be amenable to experimental 
observation. 
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The computed quartic force field of BeH 2 is available for download, in a variety of 
coordinate systems, at http://theochem.weizmann.ac.il/web/papers/BeH2.html 
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Table 2. Equilibrium bond distance (A) and harmonic frequencies (cm 1 ) for BeH : 
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